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Abstract 


We  propose  two  nonparametric  tests  of  the  «issumption  that  imperfectly  re- 
p£iired  systems  are  minimally  repaired  in  the  models  of  Brown  and  Proschan 
(1983)  and  Block,  Borges,  and  Savits  (BBS)  (1985).  The  large  samiple  theory 
for  these  tests  is  derived  from  the  as3rmptotic  joint  distribution  of  the  survival 
function  estimator  of  Whitaker  and  Samaniego  (1989)  and  the  ordinary  em¬ 
pirical  survival  fimction  bcised  on  the  initial  failure  times  of  new,  or  perfectly 
repaired  systems.  Simulation  results  are  also  provided  for  the  null  hypothesis 
case,  and  under  the  alternatives  proposed  by  Kijima  (1989).  Models  assuming 
minimal  repair  specify  that  upon  repair,  a  fctiled  system  is  returned  to  the 
working  state,  while  the  effective  age  of  the  system  is  held  constant;  that  is, 
the  distribution  of  the  time  until  the  next  failure  of  the  repaired  system  is  the 
same  as  for  a  system  of  the  same  age  which  has  not  yet  failed.  These  models 
are  common  in  the  literature  of  operations  research  and  reliability,  and  prob¬ 
abilistic  resvdts  and  the  recently  proposed  inferential  procedures  of  Whitaker 
and  Samaniego  (1989)  and  Hollander,  Presnell,  and  Sethuraman  (1989)  de¬ 
pend  on  the  minimal  repair  assumption.  Though  tests  have  been  proposed 
for  goodness  of  fit  of  the  model  when  a  particular  form  of  the  distribution  is 
2issiimed,  we  know  of  no  previous  proposal  of  a  nonparametric  method  to  test 
this  assumption. 

KEYWORDS:  Imperfect  repair,  repairable  system,  nonhomogeneous  Poisson 
process,  martingale. 


1.  INTRODUCTION 


Models  which  allow  for  imperfect  repair  are  clearly  needed  for  reliability  anal¬ 
ysis  of  repairable  systems,  and  minimal  repair  provides  a  mathematically 
tractable  approach  to  this  modeling  problem.  Models  involving  minimal  repair 
appear  frequently  in  the  literature  of  reliability  and  operations  research  (see 
Ascher  and  Feingold  (1984)  for  references).  A  particular  exeunple  is  the  age- 
dependent  minimal  repair  model  (BBS  model)  of  Block,  Borges,  and  Savits 
(1985),  which  generalizes  the  imperfect  repair  model  of  Brown  and  Proschan 
(1983).  In  the  BBS  model,  a  system  with  survival  time  distribution  F  is  put  on 
test  at  time  zero.  Upon  failure  at  age  t,  one  of  two  types  of  repair  is  performed: 
with  probability  p{t)  a  perfect  repair  is  performed  £ind  the  system  is  returned 
to  the  “good-as-new”  state;  otherwise,  a  minimal  repadr  is  performed  and  the 
device  is  returned  to  the  working  state,  but  is  only  ais  good  as  a  working  item 
of  age  equal  to  the  age  of  the  device  at  failure.  In  the  former  case,  we  consider 
the  effective  age  of  the  system  to  be  returned  to  zero  (a  new  item  has  effective 
age  0),  while  in  the  latter  case  the  effective  age  of  the  system  is  unchanged 
from  the  effective  age  at  failure.  Thus,  if  a  minimal  repair  is  performed  on  an 
item  faiiling  at  age  i,  then  the  repaired  item  has  survival  function 


F{3  +  t) 

Fit)  ’ 


5  >  0, 


where  we  use  the  notation  F  to  indicate  the  survival  function,  1  —  JP,  of  a 
distribution  F.  The  process  is  continued  after  repair,  with  each  subsequent 
failure  being  followed  by  a  perfect  repair  with  probability  p(t),  or  a  minimal 
repair  with  probability  qit)  =  1  —pit),  where  t  is  the  effective  age  of  the  failed 
device. 

Whitaker  and  Samaniego  (1989)  and  Hollander,  Presnell,  and  Sethuraman 
(HPS)  (1989)  propose  several  inferential  procedures  for  this  model.  These 
procedmes  and  their  asymptotic  distribution  theory  depend  strongly  on  the 
assumption  that  imperfectly  repaired  systems  are  minimally  repaired,  although 
in  practice  this  a^sumpticn  is  tenuous.  Minimal  repair  is  generally  motivated 
as  an  approximation  to  the  situation  where  some  failing  component(s)  of  a 
complex  multi- component  system  is  replaced  or  repaired,  but  we  would  not 
expect  the  model  to  fit  exactly  the  repairs  of  an  actual  system.  This  is  dis¬ 
cussed  by  Bergman  (1985),  who  distinguishes  between  statistical  minimal  re¬ 
pair  (ais  defined  above)  and  physical  minimal  repeiir,  in  which  the  failed  system 
is  restored  to  its  exact  physical  condition  just  before  failure.  Some  interesting 
arguments  ageiinst  the  indiscriminate  use  of  minimal  repaur  models  are  also 
given  by  Arjas  and  Norros  (1989)  and  Natvig  (1990),  who  refer  to  the  type 
of  minimal  repair  considered  here  as  ‘black  box’  minimal  repair,  and  contrast 
this  with  models  which  take  into  account  available  information  about  the  con¬ 
dition  of  the  system.  In  spite  of  the  widespread  use  and  criticism  of  minimal 


repair,  however,  we  know  of  no  previous  proposal  of  a  method  for  testing  the 
assumption  of  minimal  repair  which  does  not  assume  a  particular  form  for  the 
underlying  survival  distribution  of  the  system. 

In  this  paper  we  propose  two  nonparametric  tests  of  the  minimal  repeiir  as¬ 
sumption  in  the  BBS  model.  These  tests  are  based  on  the  notion  of  contrasting 
two  available  estimators  of  the  system  survival  function.  We  assume  that  n 
systems  are  observed  under  the  BBS  model,  each  until  the  time  of  its  first 
perfect  repair  (or  equivalently,  a  single  system  is  observed  until  the  time  of  its 

perfect  repair).  Regardless  of  the  mode  of  repair,  a  consistent  estimator  of 
the  system  survivaJ  function  is  avciiiable,  namely  the  empirical  survival  func¬ 
tion,  Fc,  based  on  the  initial  failure  times  of  the  systems  under  observation. 
Under  the  BBS  model,  the  estimator  of  Whitaker  cind  Seimaniego  (1989),  F, 
also  provides  a  ^.onsistent  estimator  of  the  system  survival  function,  while  if 
the  minimal  repair  assumption  of  the  model  fails  to  hold,  this  estimator  might 
diverge  from  the  true  system  survival  function.  In  Section  2.  we  propose  a 
Kolmogorov-Smirnov  type  test  based  on  the  maximum  absolute  difference  be¬ 
tween  F  and  F*,  and  a  test  based  on  a  Mann- Whitney- Wilcoxon  type  statistic 
of  the  form  /  Fe  dF.  Theoretical  justification  of  these  procedures  is  given  and 
they  are  applied  to  the  data  of  Proschan  (1963).  In  Section  3.  the  results  of 
a  simulation  study  of  these  procedures  are  given.  Finally,  in  Appendix  A.  we 
give  proofs  of  the  theoretical  results  presented  in  Section  2.. 

2.  TWO  TESTS  OF  MINIMAL  REPAIR 
2.1  The  Model 

The  following  model  for  sampling  from  an  age-dependent  minimal  repair  pro¬ 
cess  is  given  in  (HPS)  and  is  repeated  here  for  convenience.  To  simplify  the 
exposition,  we  assume  that  F  is  continuous  and  that  F{t)  <  1  for  all  t  <  oo, 
although  we  generally  give  expressions  applicable  to  the  more  general  case. 

A  sequence  of  failure  ages  obtained  under  a  model  of  perennial  minimal 
repair  may  be  defined  as  follows.  Let  F  be  a  life  distribution,  and  let  {Xq  = 
0,Xi,X2,...}  be  a  record-vcilue  sequence  based  on  F;  that  is,  (Xfc)^i  is  a 
Markov  process  with  the  conditional  distribution  of  Xk  given  Xo,...,Xfc_i 
being  given  by  F{t\Xk-i)  =  F{t)/F{Xk-i),  for  t  >  Xk-i  and  fc  >  1.  Note 
that  the  counting  process 

X*(<)  =  #{fc:Xfc<<} 

is  just  a  nonhomogeneous  Poisson  process  with  mean  function  equal  to  the 
cumulative  hazard  function  of  F, 


Perfect  repair  is  now  introduced  through  the  use  of  independent  uniform  (0, 1) 
random  variables,  { I7i,  C/j, . . . },  <ilso  independent  of  (Xfc)fcLi.  Defining  Sk  = 
I{Uk  <  p(Xfc))  and  v  =  inf  { fc  :  =  1 },  where  we  take  inf0  =  oo,  we 

see  that  P{Sk  =  =  p{Xk).  Clearly,  observation 

of  {(Xi,  5i), . . . ,  (X^, is  equivalent  to  observation  of  a  system  under  the 
BBS  model  until  the  time  of  the  first  perfect  repair,  although  u  may  in  general 
be  infinite. 

Let  H  be  the  subdistribution  function  defined  by 

.ff(t)  =  P(X„<t,»/<oo). 

Then  H  has  cumulative  hazard  function 

^h(<)=  /  p(3)dA(s),  (2.1) 

Jo 

and  H  is  a,  proper  distribution  function  if  and  only  if 

f  p(s)  dA{s)  = +00 .  (2-2) 

Jo 

In  this  case,  X„  is  sdmost  surely  finite  with  distribution  fimction  H.  We  will 
assume  that  (2.2)  holds  for  the  remainder  of  the  paper. 

2.2  Theoretical  Results 

Now  suppose  that  we  observe  n  such  systems,  so  that  we  have  independent 
record  value  sequences  {  {Xjk)kLi  :  1  <  i  <  n}  from  F,  independent  uniform 
random  variables  {Ujk  and  we  observe  {(Xjfc,6,fc)  : 

1  <  k  <  I'j,  1  <  j  <  n}.  Let  F*  be  the  empirical  cumulative  distribution 
function  formed  in  the  usual  way  from  the  initial  failure  ages  of  the  system, 
Xii , . . . ,  X„i ,  let  H  be  the  empiricaJ  c.d.f.  of  Xj,^ ,  •  •  • ,  Xm/„ ,  and  let 


«■;(()  =  #{ A: 

(2.3) 

(2.4) 

(2.5) 

Note  that  N{t)  is  the  total  number  of  observed  failures  through  time  t,  and 
nll{t—)  is  the  number  of  systems  still  under  observation,  or  “at  risk,”  just 
before  t.  Let  T  be  the  first  failure  age  at  which  only  one  item  is  at  risk: 


T  =  sup{  i  :  9{t)  =  i  }  =  mm{  X,., :  f  (Afp,-)  =  i  } , 

71  71 


3 


where  the  X(i)  are  the  ordered  values  of  the  combined  observed  failure  ages  of 
the  n  BBS  processes.  A  natural  estimator  of  the  cumulative  hazard  function 
is  then 

''0  nH{s—)  x^i^KtAT  nH{X(i)—) 

where  J{s)  =  I{s  <  T),  and  we  take  §  =  0.  With  this  notation,  the  survival 
function  estimator  of  Whitaker  and  Samaniego  (1989)  is 


= n(i  -  aaw)  =  n  (i-  ■■)■ 

x^i)<tAT 


t<t 


We  define 


c«)=/; 


dF(,,) 


A  H(s-)F(s) 


and 


G  =  L/{1  +  L),  G=i/(1  +  L). 

The  following  results  are  proven  in  Appendix  A.: 

Theorem  2.1.  As  n  — »  oo, 

^  B(L)  inDIO.oo). 


(2.7) 


(2.8) 

(2.9) 

(2.10) 


where  B  is  Brownian  motion  on  [0,  oo). 


Corollary  2.2. 

^  B"((?)  inD|0,oo), 

and 

v^£(F-.F.)  B°(G)  inl>[0,oo), 

F 

where  is  the  Brownian  bridge  on  [0, 1]. 

Corollary  2.3. 

y/^{F-F^)  FxB{L)  inD[0,oo]. 


4 


mo. 


Corollary  £.4.  Let  V  =  /~  F*  dF.  Then 


where 


1 

12 


f(») 

H(,-) 


dF(s) . 


(2.11) 


The  result  of  Corollary  2.2  leads  to  a  Kolmogorov- Smirnov  type  test  of  the 
minimal  repair  assumption,  to  be  carried  out  by  referring  the  statistic 


Sr  =  sup  -  ■?.(«)! 

0<t<T  F{t) 

to  a  table  of  the  distribution  of  the  supremum  of  the  absolute  Vcilue  of  the 
Brownian  bridge  over  the  interval  [0,  (?(t)]  (Koziol  and  Byax  1975;  Hall  and 
Wellner  1980).  Corollary  2.4,  on  the  other  hand,  suggests  that  we  refer  the 
statistic 

y  =  jav 

to  a  standard  normal  distribution,  where 


=1-1  f 

^  12  4/0 


^3 

~  F  is) 


4^^dFis). 

H{s-) 


(The  consistency  of  a  is  demonstrated  in  HPS,  pp.  27-28.)  Other  tests  based 
on  the  simple  idea  of  searching  for  a  statisticedly  significant  discrepancy  be¬ 
tween  Fe  and  F  are  also  possible,  but  we  will  consider  only  these  two. 


Remark  2.5.  A  strictly  monotone  transformation  of  the  time  axis  does  not 
alter  the  value  of  V*  or  the  value  of  5r,  if  t  is  also  transformed.  But  if  p{t)  =  p 
is  a  constant  function,  then  the  transformation  Xjj^  =  A(Xjk)  reduces  us  to 
the  situation  of  observing  the  BBS  model  with  the  same  constant  p(-)  and 
with  exponential  F.  Thus,  when  p(')  is  a  constant  function,  the  distribution 
of  V*  under  the  null  hypothesis  does  not  depend  on  F.  If  u  is  fixed  and 
T  =  F~^(u),  then  the  null  distribution  of  Sr  is  also  independent  of  F.  Of 
course  the  distributions  in  both  cases  depend  on  p  and  in  the  latter  case  on  u. 


2.3  An  Example 

As  an  example,  we  have  applied  the  procedures  above  to  the  Boeing  eur  condi¬ 
tioner  data  of  Proschan  (1963).  The  data  are  reproduced  in  Table  1.  For  this 
analysis,  we  have  treated  the  intervals  between  failures  as  inter-fciilure  times 
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7907 

7908 

7909 

7910 

7911 

Plane  Number 
7912  7913  7914 

7915 

7916 

7917 

8044 

8045 

T54“ 

413 

90 

74 

55 

23 

^7^ 

50 

359 

50 

130 

487 

102 

15 

14 

10 

57 

320 

261 

51 

44 

9 

254 

493 

18 

209 

41 

58 

60 

48 

56 

87 

11 

102 

12 

5 

100 

14 

29 

37 

186 

29 

104 

7 

4 

72 

270 

283 

7 

57 

33 

100 

61 

502 

220 

120 

141 

22 

603 

35 

98 

54 

181 

65 

49 

12 

239 

14 

18 

39 

3 

12 

5 

32 

9 

14 

70 

47 

62 

142 

3 

104 

85 

67 

169 

24 

21 

246 

47 

68 

15 

2 

91 

59 

447 

56 

29 

176 

225 

77 

197 

438 

43 

134 

184 

20 

386 

182 

71 

80 

188 

230 

152 

36 

79 

59 

33 

246 

1 

79 

3 

27 

201 

84 

27 

«« 

21 

16 

88 

130 

14 

118 

44 

** 

15 

42 

106 

46 

230 

59 

153 

104 

20 

206 

5 

66 

34 

29 

26 

35 

5 

82 

5 

61 

31 

118 

326 

12 

54 

36 

34 

18 

25 

120 

31 

22 

18 

156 

11 

216 

139 

67 

310 

3 

46 

210 

57 

76 

14 

111 

97 

62 

26 

71 

39 

30 

7 

44 

11 

63 

23 

22 

23 

14 

18 

13 

34 

62 

11 

191 

14 

16 

18 

130 

90 

163 

208 

1 

24 

70 

16 

101 

52 

208 

95 

Table  1.  Intervals  Between  Failures  of  Boeing  Air  Conditioner  Systems. 


between  minimal  repairs.  Proschan  omits  any  failure  interval  immediately  fol¬ 
lowing  a  major  overhaul  (indicated  by  **  in  Table  1),  and  we  have  omitted  the 
intervals  following  such  an  overhaiil  from  our  analysis,  since  it  is  impossible 
to  determine  the  age  of  the  unit  after  the  overhaul.  Thus  the  values  below 
the  **’s  in  Table  1  are  not  used.  We  have  treated  the  age  at  which  a  major 
overhaul  occurs  as  the  time  of  the  first  perfect  repair  for  that  airplane.  This 
affects  planes  7908,  7909,  7910,  and  7911.  For  purposes  of  this  exeumple,  we 
treat  the  last  observed  failure  ages  of  the  remaining  planes  as  the  times  of 
their  first  perfect  repair. 

We  have  arbitrarily  decided  to  compute  the  S  statistic  over  the  interval 
from  0  to  500  hours.  The  value  of  S  in  this  case  is  0.7705,  and  G(500)  = 
.9902.  This  value  of  5  is  less  than  the  50th  percentile  of  its  asymptotic  null 
distribution,  so  that  the  test  yields  no  evidence  against  the  minimal  repair 
assumption.  For  the  Wilcoxon-like  test  statistic,  V  =  0.4984,  dy  =  0.1753, 
and  V*  =  —0.03323,  again  yielding  no  evidence  against  the  minimal  repair 
assumption.  In  ’T'igure  1,  we  have  plotted  F  and  F*  on  the  interval  [0, 500].  In 
agreement  with  the  results  of  the  test  procedures,  there  is  little  visual  evidence 
to  suggest  a  meaningful  discrepancy  between  the  estimators.  As  we  shall  see 
in  Section  3.,  however,  these  tests  appear  to  have  very  low  power  for  ?  sample 
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size  as  small  as  13,  so  that  a  lairger  sample  size  is  needed  to  effectively  test  the 
minimal  repair  assumption  in  this  example. 

3.  SIMULATION  RESULTS 

We  have  carried  out  simulation  studies  of  both  the  size  and  power  of  the 
tests  proposed.  Computations  were  carried  out  on  the  University  of  Florida 
Statistics  Department’s  system  of  DEC  workstations.  Simulation  programs 
were  written  in  FORTRAN.  The  uniform  generator  used  as  a  basis  for  all 
other  random  numbers  was  the  imiversal  random  n\imber  generator  described 
in  Marsaglia,  Tsang,  and  Zaman  (1990).  Exponential  random  variables  were 
generated  using  the  fxmction  REXP  given  by  Marsaglia  and  Tsang  (1984).  The 
memoryless  property  of  the  exponential  distribution  meikes  it  eeisy  to  gener¬ 
ate  record  values  by  simply  adding  an  independent  exponential  to  the  current 
record  value.  For  the  general  gamma(a)  distribution,  the  initial  record  val¬ 
ues  were  generated  in  the  obvious  way:  generate  gamma  v£U’iates  until  the 
current  record  is  exceeded.  These  gamma  random  variables  were  generated 
using  the  squeeze  method  of  Marsaglia  (1977).  Once  this  sequence  produced 
a  value  greater  th2tn  a  (arbitrary,  but  it  is  clear  that  at  some  point  the  start¬ 
ing  algorithm  must  be  abandoned),  an  algorithm  described  by  Marsaglia  and 
Tsang  (1984)  f  >r  generating  random  variables  from  the  tail  of  a  distribution 
was  employed. 

The  sample  sizes  examined  were  10,  20,  30,  50,  100,  and  200.  For  each 
sample  generated,  both  S  and  V*  were  calculated.  Since  the  interval  [0,r] 
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V* 


Figure  2.  Size  of  S  and  V*  with  nominal  size  0.10,  with  p  held  constant  at 
0.50  ( — ),  0.25  (-),  and  0.10  (solid).  Estimates  based  on  30,000  iterations. 

on  which  the  statistic  is  computed  may  affect  the  behavior  of  the  S  test,  we 
collected  results  of  computing  the  statistic  to  the  90**^,  95*^,  and  99*^  percentile 
of  the  underlying  distribution  By  a  slight  abuse  of  our  earlier  notation,  we 
will  denote  these  statistics  by  5.9o,  etc.  The  tables  of  Koziol  eind  Byar  (1975) 
were  used  to  find  the  appropriate  value  of  Aa(if(T)),  after  first  rounding  K(t) 
to  the  nearest  tenth.  The  BBS  process  was  simulated  with  constcint  p  taking 
values  .10,  .25,  and  .50. 

The  estimated  sizes  of  the  tests  based  on  5.95  and  V*  with  nominal  size 
a  =  0.10  are  presented  graphically  in  Figure  2.  (The  results  for  S.90  and  5.99 
are  similar  to  those  for  5.95.)  These  estimates  are  based  on  30,000  iterations  of 
the  simulation.  Recall  that  under  the  null  hypothesis  with  constant  p,  the  null 
distributions  of  the  statistics  considered  here  do  not  depend  on  the  underlying 
F.  The  light  dotted  lines  at  a  =  0.10  ±  1.96  x  O.3/-v/30,0OO  are  provided 
for  reference.  Both  the  5,95  and  the  V*  i  ?sts  are  conservative  at  the  sample 
sizes  examined,  with  the  discrepancy  between  the  nominal  and  the  true  sizes 
decreasing  monotonically  with  the  sample  size.  The  value  of  p  also  has  an  effect 
on  the  size,  particularly  in  the  case  of  V*,  with  larger  p  yielding  a  smaller  size 
for  the  test.  Since  the  tests  are  at  least  not  anticonservative,  we  continue  by 
examining  their  power. 

To  study  the  power  of  the  test  procedures,  we  consider  two  alternatives  to 
minimal  repair.  Models  1  and  2  of  Kijima  (1989).  Let  Xk  represent  the  actual 
age  (as  measured  by  the  clock)  and  let  Yk  represent  the  effective  age  (Kijima 
calls  this  the  virtual  age)  of  the  system  at  the  fc^^fjulure.  Thus  the  conditional 
survival  function  of  AXfc  =  Xk  —  Xk-i  given  {  Xi, . . . ,  Xk-iyVi, . . . ,  Yk-i  }  is 
F(s|  Ifc-i)  =  P(s+yfc_i)/F(yfc_i),  for  s  >  0.  In  Kijima’s  Model  1,  the  effective 
age  of  the  system  after  the  fc^^'repair  is 


n  =  n-i  +  AkAXk , 


where  Ak  is  a  random  variable  taking  values  in  the  unit  interval  independently 
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of  Vi, ,  Vfc-i  and  Xi, . . . ,  Xfc.  The  idea  here  is  that  the  repair  affects  only 
aging  of  the  system  accumulated  since  the  last  repair.  In  Kijima’s  Model  2  on 
the  other  hand,  the  effective  age  of  the  system  after  the  fc*^repair  is  given  by 


n  =  ^*(n-i  +  Ax*), 


where  A}t  is  as  in  Model  1. 

In  our  simulation  study,  the  AkS  in  these  models  were  tciken  to  be  const2uit 
(vauriously,  0.75,  0.50,  0.25,  and  0.10).  The  underlying  distributions  considered 
were  gaimma  with  shape  parameters  2,  4,  and  6.  In  lieu  of  tables  of  the 
simulation  results,  we  provide  Figures  3-8.  Again  we  provide  results  only  for 
5.95,  as  the  results  for  S.90  and  S.99  were  similar.  These  figures  support  the 
following  conclusions: 

•  The  power  of  both  tests  increases  as  Ak  decreases.  This  is  to  be  expected 
since  values  of  Ak  close  to  one  in  either  model  correspond  to  models  which 
are  “close”  to  the  minimal  repair  model. 

•  With  edl  other  factors  being  constant,  the  power  of  both  tests  tends  to 
be  higher  for  Model  2  than  for  Model  1.  This  is  not  surprising,  since  for 
the  same  value  of  Ak,  repair  under  Model  2  has  a  greater  effect  on  the 
effective  age  than  repadr  under  Model  1. 

•  The  power  of  both  tests  increases  as  p  decreases.  Again,  this  is  expected, 
since  smaller  values  of  p  correspond  to  the  observation  of  more  repairs 
and  hence  more  information  about  the  mode  of  repeiir. 

•  The  power  of  both  tests  increases  as  the  gamma  shape  parameter  a 
increases  away  from  1.  This  again  is  not  surprising,  since  values  of  a  close 
to  1  correspond  to  an  underlying  distribution  closer  to  the  exponentied, 
for  which  both  Models  1  and  2  are  equivalent  to  the  minimal  repair 
model. 

•  The  test  based  on  V*  tends  to  have  greater  power  for  large  sample  sizes 
(>  30),  while  the  test  based  on  5,9$  has  greater  power  for  small  sample 
sizes  (<  30).  We  have  no  obvious  explanation  for  this  behavior. 

Since  we  know  of  no  competing  procedures  which  do  not  assume  a  partic- 
uletf  form  for  F,  it  is  difficult  to  make  an  overall  qualitative  judgement  of  the 
power  of  these  teste.  In  either  of  the  models  considered,  a  procedure  based  on 
estimation  of  a  parameter  of  the  distribution  of  Ak  (such  as  the  mean  when  the 
Afc’s  are  random  and  i.i.d.)  might  have  better  power,  but  might  not  perform 
as  well  against  other  types  of  edternatives  as  our  more  general  approach. 
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V* 


^S5 


Figure  3.  Estimated  power  (5000  simulations)  of  5.95  and  V*  against  Kijima’s 
Model  1  with  F  =  Gamma(2),  and  An  =  0.75(o),0.50(O),  0.25(0),  0.10(A). 


S.,5  V* 


Figure  4.  Estimated  power  (5000  simulations)  of  S,gs  and  V*  against  Kijima’s 
Model  1  with  F  =  tiamma(4),  and  An  =  0.75(o),0.50(n),  0.25(0),  0.10(A). 
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V* 


102030  SO  100  200  102030  SO  100  200 


Figure  5.  Estimated  power  (5000  simulations)  of  S,qs  and  V*  agednst  Kijima’s 
Model  1  with  F  —  Gamma(6),  and  An  =  0.75(o),  0.50(D),  0.25(0),  0.10(A). 


102030  so  100  200  102030  SO  100  200 


Figure  6.  Estimated  power  (5000  simulations)  of  Sss  and  V*  against  Kijima’s 
Model  2  with  F  =  Gamma(2),  and  An  =  0.75(o),  0.50(D),  0.25(0),  0.10(A). 
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V* 


^S5 


Figure  7.  Estimated  power  (5000  simulations)  of  5.95  and  V*  against  Kijima’s 
Model  2  with  F  =  Gamma(4),  and  An  =  0.75(o),0.50(n),  0.25(0),  0.10(A). 


Figure  8.  Estimated  power  (5000  simulations)  of  8.9$  and  V*  against  Kijima’s 
Model  2  with  F  =  Gamma(6),  and  An  =  0.75(o),  0.50(D),  0.25(0),  0.10(A). 
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A.  PROOFS 


In  this  section,  we  follow  the  usual  convention  of  taking  all  (7-fields  to  contain 
the  P-null  sets  of  T.,  where  (Q,  P)  is  our  basic  probability  space.  Most  of  the 
results  on  stochastic  integration  which  we  need  here  are  conveniently  collected, 
with  references,  in  Appendix  B  of  Shorack  and  Wellner  (1986):  in  particular, 
we  will  generally  use  Theorem  B.3.1  without  further  citation.  In  this  regard, 
we  note  here  that  for  all  the  local  martingales  below,  the  localizing  sequence 
of  of  stopping  times  may  be  taken  to  be  any  sequence  of  constants  diverging 
to  positive  infinity  (Hollander  et  al.  1989,  Remark  2.2).  Also,  integrands  are 
readily  verfied  to  be  predictable,  and  generally  are  bounded  on  finite  intervals 
by  some  constant,  so  that,  e.g.,  condition  (b)  of  this  theorem  is  easily  verified. 
Additioned  results  which  we  will  need  are  the  related  Proposition  18.13  of 
Metivier  (1982)  and  Lemmas  A.l  and  A.2  of  Doss  and  Chiang  (1991).  Since 
the  localizing  sequences  for  the  processes  considered  may  be  taken  to  be  a 
common  sequence  of  constants,  we  c«in  use  the  latter  two  lemmas  (which  eire 
stated  for  martingales)  without  modification. 

A.l  The  Basic  Martingale  Dependence  Structure 

Returning  to  the  basic  model  of  Section  2.,  let  Nj  be  defined  as  in  (2.3)  and 
let 

=  a({N;(3) :  s  <  t}U  :  fc  >  1}). 

Then 

M;(t)  =  «;(()  -  A(<) 

is  a  locally  square-integrable  martingale  with  respect  to  with  predictable 
quadratic  variation 

(M;)(t)=  /‘(l-AA(5))dA(s), 

V  0 

and  with  Nj  as  in  (2.4), 

M,(t)  =  Ar,(<)-A(tAX,v,)=  f  I{Xi^,>s)dM*{s) 

Jo 

is  a  locally  square-integrable  .7^'^^-martingale  with 

{Mi){t)  =  /‘/(X,v,  >  s)(l  -  AA(s))(fA(5) 

J  0 

(Hollander  et  al.  1989).  Similarly,  we  note  that  I{Xj\  <t)  =  Nj{t  AXji)^  and 
that  (note  that  both  I{Xji  >  i)  and  I{Xj„-  >  t)  are  predictable), 

M;(t)  =  I{Xjr  <t)-  A(t  A  Xj^)  =  /(X,i  >  s)  dM;{s) 
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is  a  square-integrabie  ^^^-martingale  with 

(M;)(t)  =  f  I(Xir  >  s)(l  -  AA(s))iA(s). 

JO 

Fin£illy,  by  Proposition  18.13  of  Metivier  (1982), 

=  /‘/(A-^.,>s)/(Jfy.>s)d(A/;)(s) 

JO 

=  r  liXjx  >  s){l  -  AA(s))  dA{s).  (A.l) 

Jo 

Now  let 

i=i 

Then,  noting  that  nPjJt)  =  Y!i=\  KXjx  <  <)>  =  Ei=i  >  t), 

and  nH(t—)  =  —  ^)>  recalling  the  definition  of  N  in  (2.5), 

it  follows  from  the  independence  of  the  (Doss  and  Chiang  1991,  Lemmzis 
A.l  and  A. 2)  that 

Afe(<)  =  E  M;(<)  =  nF,(t)  ~nf  f,{s-)  dfx{s)  (A.2) 

i=i 

and 

M{t)  =  E  -  "  /*  ^(^-)  (A.3) 

i=i 

are  locally  square-integrabie  .7^t*niartingales  with 

(M,)(t)  =  =  n  /‘f.(s-)(l  -  AA(s))dA(s) 

,=l 

and 

(Af)(t)  =  X:(M,)(()  =  n  /'#(s-)(l  -  AA(s))  JA(s). 

i=i  •'« 

Since  M,-  and  Mj,  are  independent  when  j  ^  f,  we  may  agedn  apply  Lemma 
A.2  of  Doss  cind  Chiang  (1991)  to  get  {Mj^My)  =  0.  Thus,  by  (A.l)  and  the 
bilineaurity  of  the  (*,  •)-operator 


3=X  j’=X  j=X 

=  Er/(X,i>s)(l-AA(s))dA(5) 

j=X  •'0 

=  n  /‘  Fe(s-)(1  -  AA(3))  dA{s). 

Jo 


(A.4) 
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These  resTilts  provide  us  with  the  basic  martingale  structure  needed  to 
prove  the  results  of  Section  2..  The  first  step  is  to  establish  a  martingale 
representation  for  (F  —  F^)fF. 


A.2  Martingale  Representations 

Recalling  (2.6)  and  (2.7),  we  define 

2md  note  that 


$<t 


For  brevity,  let  Z  =  y/n{^F  —  F)IFy  =  ■)/n{Ft  —  F)/ F,  and  Z  —  Z  —  Ze  — 
y/n(F  —  Fe)/F.  Using  Duhammel’s  equation  (Gill  and  Johansen  1990;  or  see 
Shorack  and  Wellner  1986,  Lemma  7.2.1)  and  (A.2),  we  find  that 


A 

z.{t)  = 

1  1 


=  —  / 
■v/n  Jo 


dM,{s), 


y/n  Jo  F{s) 

and  similarly  (Hollander  et  al.  1989), 


(A.5) 


Z{t)  = 

=  r  <^M(.0- 

y/nJo  H{3-)F{s) 


(A.6) 


In  both  (A.5)  and  (A.6),  the  integrand  on  the  right-hand_side  is  predictable 
and  bounded  on  any  finite  interval  (0,t]  by  a  constand  (1/F(t)).  Thus  Z,  Z*, 
and  Z  are  locally  square-integrable  martingales  with 

(z.)w  = 

nJo  F  (s) 


and 


=  /‘f  f. 

Jo  ^  f\s-)F{3) 

{Z)(t)  =  -  /*  )  d(M)( 

nJo  \H{s-)Fis)) 

=  /‘ZV 

Jo  H(.- 


f«  f\^-)  _  dF{s)^ 
#(3_)  f\s-)F{s) 


15 


Also,  applying  Proposition  18.13  of  Metivier  (1982)  to  (A.4),  (A. 5),  and  (A.6), 
we  see  that 


ft  Fis-)hi3-)  dFjs) 

Jo  H{s-)  F\s-)F{s)  ■ 


These  results  combine  to  give 


(Z)(()  =  (Z)(<)  +  (Z40)-2(z.z.)W 


dF(s) 

F\s-)F(s)  ■ 


(A.7) 


The  martingale  representation  for  Z  will  be  used  below  with  Rebolledo’s 
Martingade  Central  Limit  Theorem  to  prove  Theorem  2.1.  In  particuleu-,  (A.7) 
will  be  used  to  determine  the  variance-covariance  function  of  the  of  the  limiting 
Gaussian  process  found  in  the  theorem. 


A. 3  Proof  of  Theorem  2.1 

By  the  Glivenko-Cantelli  Theorem,  F*  a^d  H  are  uniformly  strongly  consistei^ 
for  F  and  H  respectively.  Similarly,  F  is  uniformly  (strongly)  consistent  for  F 
(Whitaker  and  Samaniego  1989;  Hollander  et  al.  1989).  Thus  it  follows  from 
(A.7)  that  for  t  <  oo. 

The  result  of  Theorem  2.1  will  follow  upon  verification  of  Rebolledo’s  strong 
ARJ(2)  condition  for  Z.  For  an  arbitrary  local  martingade  m  and  e  >  0,  let 

(T‘[m](t)  =  53  |Am(s)|V(lAm(s)|  >  e). 


«<( 


and  let  a*[m]  be  the  compensator  of  <r*[m].  Then  we  must  show  that  cr*[Z](<) 
converges  in  probability  to  0  for  each  t  <  oo.  For  this,  we  first  note  that  by 
(A.2),  (A.3),  (A.5),  and  (A.6), 


and 


AZ(s)  =  ~  AM(a)  =  Aj\r(s) , 

^F{s)H{s-)  V^F{s)H{s-) 
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Thus  the  jumps  of  and  Z  are  postitive,  and  since  |AZ(s)|  =  |AZ(s)  — 
AZe(s)|,  we  have 

|A2(»)|^7(|a2(^)|  >  £) 

<  I AZ(i )|V(|AZ(i)|  >  £)  +  |AZ,(»)l'/(lAZ.(i)|  >  6) . 


Prom  this  inequality  we  see  that 

a'[Zl(()  +  <r*|Z.l(l)-a‘[Z|(f) 

=  £(|AZ(»)|»/(1AZ(.)|  >  £)  +  |AZ.(.)|=/(|AZ,(.)|  >  €) 

-|AZ(3)|»/(|AZ(»)|  >  £)) 

is  an  increasing,  nonnegative  process.  This  implies  that  5’‘[Z],  the  compensator 
of  ^'[Z],  is  boimded  by  the  sum  of  the  compensators  of  o‘*[Z]  and  (7‘[Ze].  (To 
see  this,  note  that  more  generally,  if  A  and  B  are  increasing  processes,  cind 
i4  —  B  is  increasing,  then  A  —  B  has  compensator  equal  to  A — B,  the  difference 
of  the  compensators  of  A  and  B,  respectively.  Since  the  compensator  of  an 
increasing  process  is  nonnegative,  this  implies  that  A  >  B.)  Noting  now  that 
(Hollander  et  al.  1989) 


and  similarly. 


it  follows  easily  that  d*[Z](t)  — »  0  and  5‘*[Ze](t)  — »  0  in  probability  for  all 
t  <oo.  □ 


A. 4  Proofs  of  Corollaries 

A.4.1  Proof  of  Corollary  2.2 

The  first  part  of  Corollary  2.2  follows  directly  &om  Theorem  2.1  by  a  check 
of  covairiances.  This  is  a  standard  transformation,  sometimes  ceilled  Doob’s 
tramsformation.  The  second  part  of  the  corollary  follows  from  the  first,  and 
the  uniform  consistency  (on  finite  intervals)  of  OfF.  This  consistency  follows 
easily  from  the  definition  of  G  in  (2.10)  and  the  consistency  of  F,  and  H. 


17 


A.4.2  Proof  of  Corollary  2.3 

Weak  convergence  in  D[0,  oo]  of  y/n{F  —  F^is  established  in  Hollander  et  al. 
(1989),  eind  of  course  convergence  of  y/n^Fc  —  F)  in  D[0,  oo]  is  well  known. 
This  is  enough  to  establish  tightness  of  the  joint  distribution  of  y/n{F  —  F)  and 
■y/nC-Fe  —  F)  in  (D[0,  oo])*.  Since  the  limiting  distributions  are  concentrated 
on  (7[0,  oo],  amd  since  addition  and  subtraction  are  continuous  at  continuous 
points  in  .D[0,  oo],  it  follows  that  ■y/n(-P’ ~  •?«)  is  tight  in  £)[0,  oo].  Theorem  2.1 
indentifies  the  finite-dimensional  distributions  of  the  limiting  process,  and  the 
corollary  follows. 

A.4.3  Proof  of  Corollary  2.4 

Corollary  2.4  can  be  proved  directly  in  the  same  fashion  as  the  proof  of  The¬ 
orem  5.1  of  Hollander  et  al.  (1989),  or  by  using  the  Wilcoxon  example  of  Gill 
(1989).  We  pursue  the  latter  method  here.  In  what  follows,  convergence  in 
distribution  indicates  weak  convergence  is  in  D[0,  oo]  with  the  usual  Skorohod 
topology. 

By  Corollary  3.1  of  Hollamder  et  al.  (1989),  y/n{F  —  F)  converges  in  dis¬ 
tribution  to  a  Gaussicin  process  with  covariance  function  F{s)F{t)C{s  A  t), 
while  it  is  well  known  that  y/n{Fe  —  F)  converges  in  distribution  to  a  Gaus¬ 
sian  process  with  covariance  function  F(s  A  t)F{s  V  t).  These  results  together 
with  jClorollary  2.3  identify  the  limiting  joint  distibution  of  y/n{F  —  F)  and 
y/n{Ft  —  F)  in  (£>[0,  oo])*  as  that  of  (Wi,  W\  -f  W2),  where  W\  =  FB\{C)  and 
=  FB2{L).  Here  C  and  L  are  as  defined  in  (2.8)  and  (2.9),  and  Bi  and  B2 
are  independent  Brownian  motion  processes. 

Now,  by  Lemma  3  aad  Theorem  3  of  Gill  (1989), 

P  Fd{W2  +  W2)  +  Pw^dF 

=  HWidF. 

Jo 

This  implies  that  the  limiting  distribution  of  y/n{V  —  \/2)  is  normal  with  mean 
0  and  variance 

a*=  r  r F{s)F{t)L{s  At) dF{s)dF{t) 

Jo  Jo 

which  simplifies  to  (2.11). 
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